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We report on a theoretical study of transport properties of two coupled 
Josephson junctions and compare two scenarios for controlling the current- 
voltage characteristics when the system is driven by an external biased DC 
current and unbiased AC current consisting of one harmonic. In the first 
scenario, only one junction is subjected to both DC and AC currents. In 
the second scenario the signal is split - one junction is subjected to the 
DC current while the other is subjected to the AC current. We study DC 
voltages across both junctions and find diversity of anomalous transport 
regimes for the first and second driving scenarios. 

PACS numbers: 05.60.-k, 74.50.+r, 85.25.Cp 05.40.-a 

1. Introduction 

In symmetric devices, transport can be generated by nonequilibrium 
forces which can break space or time symmetry. In mechanical systems 
like movement of a Brownian particle in a spatially periodic and symmetric 
potential, the directed motion can be induced by an external static load 
forces or by unbiased multi-harmonic forces. Another class of systems with 
broken spatial symmetries is related to ratchet systems studied intensively 
during last twenty years [I]. In the paper, we study a relatively simple sym- 
metric system which is constructed from the well-known physical elements: 
Josephson junctions. Their role in physics is invaluable and multifaceted, 
offering a rich spectrum of beneficial applications: from the definition of 
the voltage standard, through more practical devices as elements in high 
speed circuits [2], to the future applications in quantum computing devices 
[3]. We study two Josephson junctions coupled by an external resistance. 
The evolution of the system can manifest counterintuitive nature when we 
test its response to a constant external current: it can exhibit the negative 
resistance [4j. We can formulate the general question: how to manipulate 
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the system by external drivings to get optimal and desired transport be- 
haviour? To answer this question, we propose to manipulate the system by 
two combinations of external currents. 

The paper is organized in the following way: In Sec. [2| we define the 
model and provide all necessary definitions and notation. Next, in Sec. [3j 
we study the response of the system in the case when external DC and AC 
currents are applied to one junction only. In Sec. |4j we analyze the case when 
the DC current is applied to the first junction and the AC current is applied 
to the second junction. In Sec. |5] comparison of transport characteristics 
for two driving scenarios is presented. 



2. Model 

From a more fundamental point of view, we consider a system which 
consists of two coupled (interacting) subsystems and we want to uncover 
its transport properties induced by coupling between two subsystems. As 
a particular example of the real physical structure, we study a Josephson 
junction device which consists of a coupled pair of resistively shunted Joseph- 
son junctions characterized by the critical Josephson supercurrents {Id, Ic2), 
normal state resistances (i?i,i?2) and phases {(j)i,4>2) [L5j- A schematic cir- 
cuit representing the model is shown in Fig. [T] The system is externally 
shunted by the resistance and driven by two current sources Ii{t) and 
/2(i) acting on the first and second junctions, respectively. We also include 




Fig. 1. The system of two coupled Josephson junctions shunted by an external 
resistance R3 and driven by the external currents Ii{t) and l2{t)- 



into the model Johnson-Nyquist thermal noise sources Ci(t),^2(i) and ^3(t) 
associated with the corresponding resistances Ri , R2 and R^ according to 
the fluctuation-dissipation theorem. We assume the semiclassical and small 
junction regimes where the spatial dependence of characteristics can be ne- 
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glected and photon-assisted tunnelling phenomena do not contribute to the 
general dynamics. It is the regime in which the so-called Stewart-McCumber 
model [6] holds true. The range of validity of this model is discussed in detail 
in the review paper 0. 

The Kirchhoff current and voltage laws, and two Josephson relations 
allow for the full description of the phase dynamics of both junctions within 
assumed restrictions. The dimensional form of the equations of motion is 
presented in Ref. |8]. Therein, the dimensionless variables and parameters 
are defined, and the dimensionless form of dynamics is presented. Here, we 
recall only the dimensionless version of equations of motion for the phases 
(/>! = 4>i{t) and 4>2 = 4>2{t), namely. 



P2 



(j)i = Ii(t) - Id sm(j)i +a[l2{T) - /c2sin02] + \fD t!]i{t), 
al3[h{T) - Ic2 sin 02] + a[h{T) - hi sin(/)i] + 772(r), (2.1) 



where the dot denotes a derivative with respect to dimensionless time r, 
which is defined by the dimensional time t in the following way [9] 

r='-ft, (2,2) 

where 

" 'i?l + i?2+i?3 2 ^ ^ 

are the characteristic voltage and averaged critical supercurrent, respec- 



tively. All dimensionless currents in Eqs. (2.1) are in units of Ic- E.g., 
Id — 7- Ici/Ic- The parameters 

We assume that all resistors are at the same temperature T and that the 
noise sources are modelled by independent J-correlated zero-mean Gaussian 
white noises ^i{t) {i = 1, 2, 3), i.e., {(,i{t)^j{s)) = 6ij6{t—s) for i,j S {1, 2, 3}. 
The 5-correlated zero-mean Gaussian white noises and r]2{t) which 



appear in Eqs. (2.1) are linear combinations of the original noises £,i{t) and 
the resulting dimensionless noise strength reads D = AeksT /file- 

Here we would like to stress out that the dimensional equations of mo- 
tion for phases and (p2 are symmetrical with respect to the transformation 



Ri ^ R2- However, their dimensionless equivalents (2.1) are not symmetri- 



cal with respect to the change Ri ■h^ R2- It is because of: (i) the definition 
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of the dimensionless time (2.2) which is extracted from the equation of mo- 



tion for and, in consequence, (ii) the asymmetry of Vq in Eq. (2.3) with 
respect to Ri and R2 ■ 

Sometimes, it might be helpful to image the dynamics of two Josephson 



junctions described by Eqs. (2.1) as a motion of two interacting Brown- 
ian particles driven by external time-dependent forces. In this mechanical 
framework we have the following correspondence: xi = (pi, X2 = 4'2-, where 
Xi for i = 1,2 stands for the coordinate of the first and second particles, 
respectively. The main transport characteristic of such a mechanical system 
are the long-time averaged velocities vi = {(f)i) and V2 = (^2) of the first 
and second particles, respectively. In terms of the Josephson junction system 
it corresponds to the dimensionless long-time averaged voltages vi = {(pi) 
and V2 = {(P2) across the first and second junctions, respectively (from the 
Josephson relation, the dimensional voltage V = {h/2e)d(p/dt and there- 
fore d(p/dT = V/Vq). The junction resistance (or equivalently conductance) 
translates then into the particle mobility. The phase space of the determin- 
istic system ( |2.1[ ) is three-dimensional, namely {xi = (pi,X2 = (p2, 2:3 = oot}- 
Note that it is a minimal dimension for the system to display chaotic evo- 
lution in continuous dynamical systems which may be a key feature for 
anomalous transport to arise |10| ITT1 [T2l [T3] . Other aspects of dynamics of 
two coupled Brownian particles has been studied in literature |14]. However, 
experimental realizations of such systems would be difficult to construct. 

The considered system is characterized by four dimensionless material 
constants: {Ici,Ic2^ot, (3} and by the dimensionless temperature D. Addi- 
tionally, drivings Ii{t) and I2{t) are also characterized by some parameters. 
In order to reduce a number of parameters of the model we consider a system 
of two identical junctions, i.e., Ri = R2 and Id = Ic2 = 1- In such a case 



a/3 = 1 and Eqs. (2.1) takes the symmetric form 

<pi = /i(r) - sin</)i -I- alhir) - sm(p2] + ??i(t), 

<p2 = hir) - sin 02 + a[Ii{T) - sin^i] + Vd ??2(t). (2.5) 

The parameter a plays the role of the coupling constant between two junc- 
tions and can be changed by variation of the external resistance R^. The 
set of two differential equations is decoupled for a = which results with 
two independent subsystems. It is the case when i?3 — )• cxd. Note that when 
i?3 = 0, the parameter a = 1 and the system is coupled. 



3. The first scenario: DC and AC currents applied only to one 

junction 

The external dimensionless currents /i(t) and l2iT) can be modelled in 
a various way. In experiments with Josephson junctions, 'the most popu- 
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lar' three classes of currents have been apphed: DC currents, AC currents 
consisting one harmonic and AC biharmonic currents: 

/j(r) = /j + cos(wir) + 6j cos(r2jr + i = 1,2. (3.1) 

We start with the first scenario in which we apply the external current to 
the first junction only, namely, 

/i(r) = /i + aicos(wr), h{T) = {). (3.2) 



In this special case Eqs. (2.5) take the form: 

4)1 = /i — sin(/)i — asin02 + oi cos((x)r) + ?/i(t), (3.3a) 
(f)2 = all — sin(/)2 — asinc^i + aai cos(u;r) + \/15 t]2{t)- (3.3b) 

This case was considered in Ref. [15] in the context of indirect control of 
transport and absolute negative conductance induced by coupling between 
two junctions. Here, for the reader's convenience, we recall the main trans- 
port characteristics of the system but just before we'll do it let us clarify 
some technical issues. 

The above set of equations cannot be handled by known analytical meth- 
ods for solving ordinary differential equations. For this reason we have 
carried out extensive numerical simulations. We have used the stochastic 
version of Runge-Kutta algorithm of the 2"'^ order with the time step of 
10~^ • (27r/w). The initial phases 0i(O) and 02(0) have been randomly cho- 
sen from the interval [0, 27r]. Averaging was performed over 10^-10^ different 
realizations and over one period of the external driving 27r/ui. Numerical 
simulations have been carried out using CUDA environment on desktop 
computing processor NVIDIA GeForce GTX 285. This gave us possibility 
to speed up the numerical calculations up to few hundreds times more than 
on typical modern CPUs. Details on this very efficient method can be found 
in [16]. 

The voltage Vi = Vi{Ii),i = 1,2, is typically a nonlinear and non- 
monotonic function of the DC current Ii. In the normal transport regime 
the nonlinear resistance or the static resistance Ri = Ri{Ii) = Vi{Ii)/Ii (or 
equivalently conductance d = 1/Ri) is positive at a fixed bias Ii. When 
the system response is opposite to the external driving, i.e., when Ri < 
we reveal the anomalous transport regime with absolute negative resistance 
(ANR) im [T2] or nonlinear negative resistance (NNR) |13j. 

Now we would like to address some general comments about the long- 



time behaviour of the considered system (3.3). Let us consider the voltages 



^'i = ^^i(-^i) and V2 = V2{Ii) as functions of the DC bias. If we make the 



transformation Ii — t- —Ii to Eqs. (3.3), we note that vi{—Ii) = —vi{Ii) and 
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V2{—Ii) = —V2{Ii) (because the functions siiKpi and cos(a;r) are symmetric 
and noises r/i(r) are also symmetric). From these relations it follows that 
vi{0) = -t;i(0), V2i0) = -V2i0) and we deduce that t>i(0) = and i;2(0) = 
when the DC bias is zero, i.e., for Ii = 0. From results of Ref. [15j it 
follows that for high frequency (a; > 5), the long time averaged voltages Vi 
are negligible small. It is because very fast positive and negative changes of 
the driving cannot induce transport. If the DC current is sufficiently large, 
it is rather obvious that voltages across both junctions has the same sign 
as the DC bias and depend (almost) linearly on the DC current. For the 
DC current Ii > 0, one can identify three remarkable and distinct transport 
regimes: 

(I) vi> and V2 > 0, 

(II) vi> and V2 < 0, 

(III) vi < and V2 < 0. 

The regime (IV): fi < and V2 > has not been detected. From the analy- 
sis reported in Ref. [15] it follows that the most interesting transport effects 
can take place in the regime of small Ii. Indeed, it has been found that 
for Ii < 0.1 the absolute value of the long-time average voltage across both 
junctions takes its highest values for cj < 1. For stronger coupling, strips of 
non-zero average voltage begin to appear at progressively lower values of the 
amplitude ai of the AC driving. They are also visible for the average volt- 
age of the first junction, which means that the strips represent the regimes 
in the parameter space where both junctions operate synchronously. It is 
illustrated in Fig. |2] In this regime of parameters, we can detect several 
interesting effects: 

• the DC voltage vi is positive but the voltage V2 exhibits ANR for 
Ii — )• and NNR for larger value of Ii. 

• There are two different mechanisms generating negative resistance in 
the second junction: 

— In the case Ii = 0.025, the negative resistance is induced by 
thermal fluctuations: for D = the voltage ^2 > and when 
temperature increases V2 becomes negative. There is a restricted 
interval of temperature where the voltage V2 is negative. 

— In the case Ii = 0.0455, the negative resistance is generated by the 
deterministic dynamics because in the deterministic limit (when 
D = 0) the DC voltage V2 < 0. 
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Fig. 2. (color online). The first scenario: the long-time averaged voltages Vi for 
i = 1,2 of the first and second junction. Panel (a) illustrates the dependence on 
the DC current 1\ at the fixed temperature D = 2 • 10^^. In the panel (b) and 
(c) dependence on temperature D is depicted for 1\ — 0.025 and I\ — 0.0455, 
respectively. Other parameters read: coupling strength a ~ 0.77, amplitude a\ = 
1.775 and frequency w = 0.1875 of the AC driving, /2 = ft2 = 0. 



The anomalous transport effects like absolute negative resistance cannot 
occur in the decoupled system because in this case two decoupled and in- 
dependent equations correspond to the overdamped dynamics for which the 
long-time average v\ = has the same sign as Ii and V2 = (02) = 0. Fig. 
[3] shows how the average voltage across the junctions depends on the cou- 
pling constant a. One can note windows of a for which anomalous transport 
can be observed. 



4. The second scenario: DC applied to the first junction and AC 
applied to the second junction 

Nowadays technology allows experimentalists to apply driving to each 
of the junctions separately. In the following we would like to consider the 
scenario in which the DC current /i(t) = Ii is applied to the first junction 
and the AC current I2{t) = 02Cos(ci;r) is applied to the second junction. 



The corresponding dynamics is described by the special case of Eqs. (2.5), 
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Fig. 3. (color online). The first scenario: the long-time averaged voltages Vi for i — 
1,2 across the first and second junction. The infiuence of the coupling parameter 
a on transport properties is depicted in the panels (a) and (b) for two fixed values 
of the DC current Ii — 0.025 and Ii = 0.0455, respectively. Other parameters are: 
temperature D = 2- 10^^, amplitude Oi — 1.775 and frequency lo — 0.1875 of the 
AC driving, /2 = 02 = 0. 



namely, 

01 = Ii — sin0i — asin02 + 002 cos(cjr) + \/Z) ?7i(t), (4.1a) 

02 = a/i — sin02 — asin0i + 02 cos(cjr) + \/Z? ?72(''")- (4- lb) 

This scenario leads to transport characteristics which in general are different 
than in the first scenario. In particular, for the DC current Ii > 0, one can 
identify only two transport regimes where: 

(I) wi > and V2 > 0, 

(II) f 1 < and V2 > 0. 

The regimes {vi > and V2 < 0} and {vi < and V2 < 0} have not been 
detected. This means that for this type of driving the transport properties 
are a little bit modest. 

The regime (II) seems to be more interesting. In Fig. [Ij we present 
the cur rent- volt age characteristics in this regime. The unique feature is the 
emergence of the absolute negative resistance and the interval of Ii where the 
averaged voltage across the first junction vi is negative. This is to be con- 
trasted with the voltage V2 which assumes only positive values. The regime 
of the nonlinear negative resistance is not found in this scenario. The most 
profound ANR effect occurs for the dc current Ii = 0.0067. For this value, in 
panel (b) of Fig. |4j we show the voltage dependence on temperature D of the 
system. A closer inspection of the panel (b) of Fig. |4] reveals a mechanism 
responsible for generating of anomalous transport. The negative resistance 
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Fig. 4. (color online). The second scenario: the long-time averaged voltages Vi for 
i = 1,2 across the first and second junction. Panel (a) illustrates the dependence 
on the DC bias Ji at fixed value of amplitude 02 = 3.0023, frequency of AC driving 
uj = 0.1292. the coupling strength a = 0.6333 and temperature D = 4.4 • 10~^. 
Panel (b) illustrates the temperature dependence for the DC current Ii = 0.0067. 
Other parameters are the same as in panel (a) . Panel (c) shows the role of coupling 
a. Other parameters are the same as in panel (b). 



is solely induced by deterministic dynamics and even at zero temperature 
D = the resistance is negative. For this chaotic-assisted mechanism, tem- 
perature plays destructive role: if temperature increases the effect disappears 
and for temperature D greater than D = 5 ■ 10~^ the averaged voltage vi is 
positive. 



5. Comparison of transport characteristics for two scenarios 

In the previous two sections we studied properties of the DC voltage 
across the first and second junctions driven by two different external cur- 
rents. We presented the most interesting regimes where anomalous transport 
(i.e. negative resistance) can occur. There are three necessary ingredients 
for the anomalous transport to observe: DC bias, AC current (the nonequi- 
librium driving) and coupling. In this section, we compare transport prop- 
erties in the same parameter domain but for two scenarios. In Fig. |5] the 
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Fig. 5. (color online). The second scenario: the long-time averaged voltage v 
across the first (panel a) and the second (panel b) junction. The panels show the 
dependence on the dc bias 1\ at fixed temperature D = 2 ■ 10~^ for two different 
types of driving, ie. applied to only one junction (ai = 1.775,02 = 0, solid, 
blue line) and applied to both of them (oi =6,02 = 1.775, dashed, red line). 
Other parameters are coupling strength a = 0.77, the frequency of the ac driving 
Lj = 0.1875. 



current- voltage curves are compared in the regime where the ANR and ANR 
is induced for the second junction in the first scenario and the DC voltage 
across the first junction is always positive, cf. Fig. [2]^a). On the other hand, 
in the second scenario, the first junction exhibits very small absolute neg- 
ative resistance while the DC voltage across the second junction is always 
positive. It means that the unbiased AC current can change the direction 
of transport (of course, the DC current can change it but it is rather trivial 
because the DC current is biased). In Fig. |6] we compare the above charac- 
teristics in dependence on temperature (the upper panels) and the coupling 
strength (the bottom panels) in the regime where anomalous transport is 
induced by thermal fluctuations. For positive values of the DC current, one 
could expect four transport regimes: 

(I) vi> and V2 > 0, 

(II) f 1 > and V2 < 0, 

(III) vi <0 and V2 < 0, 

(IV) vi <0 and V2 > 0. 

In the first scenario, the regimes (I)-(III) can occur. In the second scenario, 
the regimes (I) and (IV) can occur. From this point of view, the first scenario 
seems to be more optimal: there are three regimes. From the symmetry of 
the system it follows that the regime (IV) could be obtained by applying 
the same driving to the second junction only. 
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Fig. 6. (color online). The long-time averaged voltage v across the first (panels 
a and c) and second (panels b and d) junction. The upper panels show the de- 
pendence on temperature D at fixed value of dc bias /i = 0.025 and the coupling 
strength a — 0.77 for two different types of driving, i.e. applied to only one junction 
(ai = 1.775, a2 = 0, solid blue line) and applied to both of them (ai = 0, a2 = 1.775, 
dashed red line). The bottom panels show the role of the coupling a for tempera- 
ture D ~ 2 ■ 10^^\ The frequency of the ac driving is a; = 0.1875. 



In summary, we studied transport properties of two coupled Josephson 
junctions and compared two scenarios for controlling the current-voltage 
characteristics when the system is driven by an external biased DC current 
and unbiased AC current consisting of one harmonic. We uncovered a reach 
diversity of anomalous transport regimes for the first and second driving 
scenarios. 
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